Dipole model analysis of the newest diffractive deep inelastic scattering data 
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We analyse the newest diffractive deep inelastic scattering data from the DESY collider HERA 
with the help of dipole models. We find good agreement with the data on the diffractive structure 
functions provided the diffractive open charm contribution is taken into account. However, the 
region of large diffractive mass (small values of a parameter 0) needs some refinement with the help 
of an additional gluon radiation. 
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I. INTRODUCTION 

Diffractive deep inelastic scattering (DDIS), observed 
at the DESY collider HERA (see []J and references 
therein) is one of the most intriguing phenomcnons in 
the electron-proton (ep) collisions. Despite high virtu- 
ality of the photonic probe, the incoming proton scat- 
ters intact being separated by a rapidity gap from a 
diffractive system, which is additionally formed in the 
final state. The understanding of these processes based 
on quantum chromodynamics (QCD) is the biggest chal- 
lenge in the area of deep inelastic scattering. In this class 
of processes, large photon virtuality Q 2 , which serves as a 
hard scale, suggests tha one use perturbative QCD with 
quarks and gluons as basic quanta. On the other hand, 
softness of the proton side and formation of the rapidity 
gap touch fundamental problems concerning transition 
into the nonperturbativc domain of QCD. Thus, such im- 
portant issues like parton saturation, unitarity and even 
confinement, arc likely to be addressed in the theoretical 
description of diffractive processes. 

The most promising QCD based approach to deep in- 
elastic scattering (DIS) diffraction is formulated in terms 
of dipole models. In these models, the diffractive, color 
singlet state is systematically built from parton compo- 
nents of the light cone wave function of the virtual pho- 
ton (see [1, HI and references therein). The lowest order 
states is formed by a quark-antiquark pair (qq) while in 
higher orders more gluons g and qq pairs are present. In 
our analysis we will concentrate on two first components, 
qq and qqg, since in the configuration space they can be 
treated as simple, quark or gluon, color dipoles. Their 
interaction with the proton is described by the dipole 
scattering amplitude N(x,r,b). Here r and b are two- 
dimensional vectors of transverse separation and impact 
parameter, respectively, and x is the Bjorken variable 
which brings the energy dependence into the dipole mod- 
els. The main advantage of this approach is the obser- 
vation that the dipole scattering amplitude can be ex- 
tracted from the DIS data on fully inclusive quantities, 
like the structure functions Fi and Ft,, based on some 
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physically motivated form with a few parameters to fit 
[E & S Hi- Then, it can be used in the description of 
diffractive processes [| M, ED, E 03 ■ The form of 
N which we use in our analysis is motivated by key fea- 
tures of parton saturation in dense partonic systems. The 
most important one is a saturation scale Q s (x) [f| which 
can be extracted from the DIS data on the structure func- 
tion i<2. The QCD based motivation for the existence of 
such a scale is provided by the analysis of the high energy 
nonlinear evolution equations of Balitsky and Kovchcgov 

[3 mm Ell- 
in this analysis, we consider two important parameter- 
isations of the dipole scattering amplitude, called Golec- 
Biernat-Wuesthoff (GBW) @ and color glass condensate 
(CGC) m which parton saturation results are built 
in. We present a precise comparison of the results of the 
dipole models which use these parametcrisations with the 
newest data from HERA on the diffractive structure func- 
tions, obtained by the HI Q and ZEUS @,Hj Collabora- 
tions. We also make a comparison with new data on the 
diffractive open charm production plj . An analysis of 
exclusive diffractive processes within the dipole approach 
was performed in [14| . Previous analyses which use par- 
ton saturation results, like those in 0, [l3|, [22|, [H, [24| . 
are based on less precises diffractive data, and in conse- 
quence, they could not address important questions re- 
lated the precise comparison presented in this paper. 

The comparison we performed prompts us to discuss 
some subtle points of the dipole models, mostly re- 
lated to the qqg component, and connect them to the 
approach based on the diffractive parton distributions 
evolved with the Dokshitzcr-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equations. Within the latter approach, 
the diffractive open charm production is particularly in- 
teresting since it is sensitive to a diffractive gluon dis- 
tribution. However, the accuracy of the existing data 
on such a production does not allow one to discriminate 
between different gluon distributions considered in our 
analysis. 

The outline of this presentation is the following. In 
Sec. [H]we present basic formulae of the color dipole ap- 
proach to diffraction while in Sec. Mil we discuss the two 
parametcrisations of the dipole scattering amplitude used 
in our analysis. In Sec. IIVI we perform a comparison of 
the dipole model results on the diffractive charm produc- 
tion with the HERA data. A similar comparison for the 
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total diffractive structure functions is presented in Sec.lVl 
In Appendix we derive a formula for the diffractive gluon 
distribution from dipole models, which is important for 
the discussion of the diffractive charm production. 



II. DIFFRACTIVE STRUCTURE FUNCTIONS 
IN DIPOLE MODELS 

In the dipole approach to DDIS, the diffractive struc- 
ture function F 2 is a sum of components corresponding 
to different diffractive final states produced by a trans- 
versely (T) and longitudinally (L) polarised virtual pho- 
ton [25[. We consider a two component diffractive final 
state which is built from a qq pair from a transverse and 
longitudinal photon anda qqg system from a transverse 
photon, see Fig.[TJ Thus, the structure function is given 
as a sum 



F, 



(qq) 



(1) 



where the kinematic variables depend on diffractive mass 
M and center-of-mass energy of the j*p system W 
through 



M 2 + Q 2 
W 2 + Q 2 ' 







Q 2 +M 2 



(2) 



while the standard Bjorken variable x = xp[3. The de- 
pendence of F 2 on the momentum transfer t = (p — p') 2 
is integrated out. The qq components from transversely 
and longitudinally polarised photons are given by 
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where / denotes quark flavours, m/ is quark mass and 
the diffractive slope Bd in the denominator results from 
the t- integration of the structure functions, assuming an 
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FIG. 1: The qq i qqg components of the diffractive structure 
function F 2 D . 



exponential form for this dependence. From HERA data, 



B d = 6 GeV- 



Thc variables 



z f = 4m}/M 2 ), Q) = z(l-z)Q 2 +m} (5) 

and the functions <f>i take the following form for i = 0, 1 

drrKi(Qfr) Ji(kfr) a(xp,r) (6) 



m 2 f is the quark transverse 



where kf = ^/z(l — z)M 2 — 

momentum while Ki and Ji are the Bessel functions. The 
lower integration limit Zf in Eqs. §3§ and Q corresponds 
to a minimal value of z for which the diffractive state with 
mass M can be produced. In such a case kf = 0. At the 
threshold for the massive quark production M 2 = 4m 2 

and Zf = 1/2, leading to = 0. For massless quarks 

The quantity a{xp, r) in Eq. © is called a dipole cross 
section and described the interaction of the qq dipole with 
the proton. It brings the energy dependence into the 
structure function formulae and is related to the imagi- 
nary part of the dipole scattering amplitude, N(xp, r, b), 
by the integral over the impact parameter 



a(xp,r) = 2 J d 2 bN(x P ,r,b) 



(7) 



Notice that for DDIS the Bjorken variable x is substi- 
tuted by xp = x/j3. For (3 ~ 1 this substitution is 
sublcading from the point of view leading logarithms of 
energy W which appear in the QCD computation of this 
amplitude. However, for large diffractive masses, (3 <C 1, 
such a substitution becomes phenomenologically impor- 
tant. 

The qqg diffractive component from transverse pho- 
tons, computed for massless quarks is given by 



xpF^ qq9) 
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where the function (f>2 takes to form 



(l-z)Q 2 \ x2 
k 2 



(8) 



dr r Ki 



1 - z 



kr J 2 (kr)a(x P ,r) (9) 



with K 2 and J 2 are the Bessel functions. In papers [9U26j|. 
formula ([5]) was computed with two gluons exchanged 
between the diffractive system and the proton. Then, 
the two gluon exchange interaction was substituted by 
the dipole cross section a = a(xp,r) for the qq dipole 
interaction with the proton. For example, for the GBW 
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paramctcrisation of the dipole cross section [26 
we discuss in the next section, is given by 
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= ct 1 - e 



2 Q 2 J4 



which fixed values of xp and Q 2 . Each component domi- 
nates in different regions of diffractive mass: dom- 
inates for M 2 ~ Q 2 ((3 ~ 1/2), F^ n) is important for 



(10) 



However, the qqg system was computed in the approxi- 
mation when parton transverse momenta fulfil the con- 
dition kr q ~ kx-q 3> kx g - Thus, in the large N c ap- 
proximation, it can be treated as a gluonic color dipole 
gg. Such a dipole interacts with the relative color factor 
Ca/Cf with respect to the qq dipole. Therefore, the two 
gluon exchange formula should be eikonalized with this 
color factor absorbed into the exponent. For the GBW 
paramctcrisation, this leads to the following gluon dipole 
cross section in Eq. j9]) 



to 1 



-(C A /C F )r 2 Q 2 j4 



(11) 



In such a case, the color factor Ca/Cf = 9/4 (for N c — 
3) disappears from the normalisation of the scattering 
amplitude and we have to rescale the structure function 
in the following way 



M19) 
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{Ca/Cf? 
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(12) 



By the comparison with HERA data, we will show in the 
next section that the latter possibility is more appropri- 
ate for the data description. 

We summarise our considerations in Fig. [2l which 
shows three components of F® as a function of f3 for 



M 2 < 



09 



1 



and F^ qq9) wins for large diffractive 



mass, M 2 > Q 2 (J3 < 1). 



III. DIPOLE SCATTERING AMPLITUDE 

We are going to compare the presented dipole descrip- 
tion of the diffractive structure functions with the newest 
HERA data. For this purpose, we consider two param- 
eterisations of the dipole cross section which are based 
on the idea of parton saturation in dense gluon systems. 
The first one is the GBW paramctcrisation with heavy 
quarks [|| which has played an inspirational role in stud- 
ies of parton saturation in the recent ten years. The 
second one is the CGC parameterisation 0, [H[ which 
somehow summarizes the studies within the Color Glass 
Condensate (2?| approach to parton saturation. Quite 
surprisingly, these two parameterisations give very sim- 
ilar results for the diffractive structure functions. The 
main reason is the same normalisation of the dipole cross 
section, Co- The origin of the same numerical value, how- 
ever, is different. For the GBW parameterisation cto is 
fitted to the data for F2 while for the CGC parameter- 
isation it is computed from a diffractive slope Bjj, see 
Eq. ([181). 

The two considered parameterisations, specified below, 
describe very well the inclusive DIS data on the structure 
function F2 . Their use for the DDIS description is a very 
important test of the universality of the dipole approach 
to DIS diffraction. 
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FIG. 2: The three components in Eq. (fl| as a function of j3 
in the massless quark limit with the GBW parameterisation 
of the dipole cross section. The qqgT component is without 
the color factor modification. 



(1) it The GBW parameterisation with heavy quarks 
has the following form of the qq dipole cross section 

1 

a(x P ,r) = ao(l-cxp(-r 2 Q2/ 4 )) (13) 

where <To = 29 vcib, and the saturation scale is given 
by 

Q 2 S = (x P /x y X GeV 2 (14) 

with xq = 4 • 1CT 5 and A = 0.288. The dipole 
scattering amplitude in such a case reads 

N(x P , r, b) = 8(b -b)(l- cxp(-r 2 Q2/ 4 ) (15) 

where 27t6q = ctq • This form corresponds to a model 
of the proton with a sharp edge. 

(2) The CGC parameterisation with heavy quarks of 
the quark dipole scattering amplitude is given by 

I III El 

N(x P ,r,b) =S{b)N(x P ,r). (16) 
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FIG. 3: The dipole cross section as a function of r for 
x = 1CT 2 . . . 10" 6 (from right to left) and for the GBW (con- 
tinuous lines) and CGC (dashed lines) parameterisations. 



where the form factor S(h) = cxp(— b 2 / (2Bd)) with 
the diffractive slope from HERA, B d = 6 GeV~ 2 . 
Thus, the dipole cross section ([7]) is given by the 
formula 



<t(xjp, r) = 4irB d N(x P ,r) 



(17) 



We see that the asymptotic value of a for r — * oo is 
the same as for the GBW paramctcrisation, if the 
diffractive slope measured at HERA is substituted, 

cr = ^B d = 29 mb . (18) 

In addition, 

N{x P ,r)= (19) 
N o( L . 

for rQ s > 2 



27 a 2 1n 2 (rQ a /2) 

e -^i-(-jp) for rQ s < 2 



1 _ e -4aln^(/3rQ a ) 



where the saturation scale Q s has now the following 
parameters: A = 0.22 and xo = 1.63 • 10~ 5 . The 
parameters a = 0.615 and (3 = 1.006 are chosen 
such that N and its first derivative are continues at 
the point r where N(r) = Nq — 0.7. The remaining 
parameters are given by k = 9.9 and 7 C = 0.7376. 

Both parameterisations provide the energy dependence 
of the diffractive structure function through the variable 
xp. This dependence is determined from fits of the dipole 
model formula for F2 into the data from HERA for the 
Bjorken variable x < 0.01. In the case of DDIS, x is 
substituted by xp. 



IV. DIFFRACTIVE CHARM QUARK 
PRODUCTION 



In the diffractive scattering heavy quarks arc produced 
in quark-antiquark pairs, cc and bb for charm and bottom, 
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FIG. 4: Top: the ccT and ccL components of from the 
dipole model with the GBW parameterisation together with 
the ccX contribution from the collinear factorisation approach 
(|23[) with the diffractive gluon distribution (IA.4|I . Bottom: 
the ccX component in a different scale against the massless 
qqT, qqL and qqg components. 



respectively Such pairs can be produced provided that 
the diffractive mass of is above the quark pair production 
threshold 



M 



1 > Am 2 , 



(20) 



In the lowest order the diffractive state consist only the 
cc or bb pair. In the forthcoming we consider only charm 
production since bottom production is negligible. The 
corresponding contributions to are given by Eqs. (J3j) 
and (|4|) with one flavour component. For example, for 
charm production from transverse photons we have 
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FIG. 5: A comparison of the collinear factorisation predictions with the GBW and CGC gluon distributions (solid lines) with 
the HERA data on the open diffractive charm production. The dashed lines are computed with the gluon distribution obtained 
in the DGLAP fit [231 to the HI data on the diffractive structure functions. 



where m c and e c are charm quark mass and electric 
charge, respectively. The minimal value of diffractive 
mass equals M^ in = 4m;?, thus the maximal value of 
(3 is given by 

Q 2 



Q 2 + Am 2 



In such a case, z c = 1/2 in Eq. (|21j) and Ff^ = for 
(3 > fimax- This is shown in Fig. [4] (top) for the cc 
diffractive states from transverse (ccT) and longitudinal 
(ccL) photons. By the comparison with the correspond- 
ing curves for three massless quarks (qqT 7 qqL), shown 
in Fig. [4] (bottom), we see that the exclusive diffrac- 
tive charm production contributes only 1 /30 to the total 
structure function F® . Thus it can practically be ne- 
glected. 

The next component is the ccg diffractive state. Un- 
fortunately, formula ([5]) for the qqg production is only 
known in the massless quark case and cannot be used for 
heavy quarks. Thus, we have to resort to the collinear 
factorisation formula, given by Eq. (|23p . in which the 
charm-anticharm pair is produced via the photon-gluon 
fusion: j*g — > cc [23|. If such an approach is applied 
to diffractive scattering, gluon is a "constituent of a 
pomeron". The diffractive state consists of additional 
particles X (called "pomeron remnant") in addition to 
the heavy quark pair, which is well separated in rapid- 
ity from the scattered proton. The collinear factorisa- 
tion formula for the charm contribution to the diffractive 
structure functions is taken from the fully inclusive case 
[2^ | in which the standard gluon distribution is replaced 
by the diffractive gluon distribution g D : 



w D(ccX) 

x P t 2L 



2(3 e 



x P g 



2 

2tt 
°(x P , 



1 dz ([3 
— G 2 ,l -, 

a/3 z \ z 



Q 2 



(23) 



where a = 1 + 4m 2 /Q 2 and the factorisation scale fi 2 = 
4m 2 with the charm quark mass m c = 1.4 GeV. The 
leading order coefficient functions are given by 

C 2 (z, r) = \ {z 2 + (1 - zf + 4z(l - 3z)r - 8z 2 r 2 } 

, 1 + a n , 
x In + \a{ - 1 + 8z(l - z) 

)r) (24) 



1 -a 
4z(l - z 



C L (z,r) 



-4z 2 r In 1 + al - a + 2az(l - z) 



(25) 



where r = m 2 /Q 2 and a = yl — 47'z/(l — z). The lower 
integration limit in Eq. (|23p results from the condition 
for the heavy quark production in the fusion: 7*5 — > cc, 

{zxpp + q) 2 > Am 2 c ( 26 ) 

where we assume that gluon carries a fraction z of the 
pomeron momentum xpp. 

The ccX contribution given by Eq. (|23|) is shown in 
Fig. [4] as the solid lines. As seen in the top figure, this 
component becomes significant for /3 < 0.1 . By a com- 
parison with the massless quark contributions (the bot- 
tom figure) we see that diffractive charm production con- 
tributes up to 30% to the diffractive structure function 
F2 for small values of (3. The presented results were ob- 
tained assuming the diffractive gluon distribution which 
results from the dipolc models, given by Eq. (|A.4[) in 
Appendix, with the GBW parameterisation of the dipole 
cross section with the color factor modification (|A.6[) . 
The CGC parameterisation gives a similar result. 

In Fig. [5] we show the collinear factorisation predictions 
for the diffractive charm production confronted with the 
new HERA data [2l| on the charm component of the 
reduced cross section: 



nice) 
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FIG. 6: The fractional charm contribution, /£> c given by Eq. (|28[) , is shown as a function of j3, for two values of xp = 0.004 
and 0.018. The solid lines are computed for the ccX contribution with the GBW and CGC diffractive gluon distributions while 
the dashed lines are found for the diffractive gluon distribution obtained in the DGLAP fit [2|| to the HI Collaboration data. 



The solid curves, which arc barley distinguishable, cor- 
respond to the result with the GBW and CGC param- 
eterisations of the diffractive gluon distributions. The 
dashed lines are computed for the gluon distribution from 
a fit to the HI data [28|] based on the DGLAP equations. 
The present accuracy of the charm data does not allow 
to discriminate between these two approaches although 
the data seem to prefer the gluon distribution from the 
DGLAP fit which is much more concentrated in the large 
z-rcgion as compared to the dipolc model gluon distribu- 
tions, see Fig. [11] in Appendix. 

The importance of diffractive charm is illustrated in 
Fig. [5] where the fractional charm contribution, 



fee „D(cc) I _D 



(28) 



to the total diffractive cross section, discussed in the next 
section, is shown as a function of (3 against the HI Col- 
laboration data [ll|. For small values of 0, the charm 
contribution equals on average approximately 20 — 30%, 
which is comparable to the charm fraction in the inclusive 
cross section for similar values of Q 2 [3fj| . 



V. COMPARISON WITH THE HERA DATA 

In Figs. [7] and [5] we show a comparison of the dipole 
model predictions with the ZEUS Collaboration data [2(| 
on the reduced cross section 



Fn 



y 



i + (i-y) 



(29) 



We included the charm contribution in the above struc- 
ture functions. The solid lines correspond to the GBW 
paramcterisation of the dipole cross section with the color 
factor modifications (fTT|) and (fT2"j) of the qqg component, 



while the dashed lines are obtained from the CGC param- 
eterisation. We see that the two sets of curves are barely 
distinguishable. This somewhat surprising results could 
be attributed to the same normalisation of the dipole 
cross section in both models, (To = 29 mb. Let us empha- 
sise again that this numerical value was obtained in two 
different ways (see Sec. IIIII for more details). The color 
factor modification of the qqg component in the GBW 
paramcterisation is necessary since the curves without 
such a modification significantly overshoot the data (by 
a factor of two or so) in the region of small (3 where the 
qqg component dominates. 

The comparison of the predictions with the data also 
reveals a very important aspect of the three component 
dipole model |T]). In the small (3 region, the curves are 
systematically below the data points, which effect may 
be attributed to the lack of higher order components in 
the diffractive state, i.e. with more than one gluon or 
qq pair. This is also seen for the HI Collaboration data 
P| shown in Fig. [9] For small values of (3 both the solid 
(GBW) and dashed (CGC) curves are below the data. It 
is also important that the charm contribution, described 
in Sec. IIV1 is added into the analysis. Without this con- 
tribution the comparison would be much worse than that 
shown here. 

This effect may be attributed to the lack of higher or- 
der components in the diffractive state, i.e. with more 
than one gluon or qq pair. They may be added in the 
DGLAP based approach to inclusive diffraction which 
sums additional partonic emissions in the diffractive state 
in the transverse momentum ordering approximation. A 
comprehensive discussion of the DGLAP based fits to the 
diffractive HERA data is presented in (28[. We only re- 
call here that in this approach the diffractive structure 
functions are twist-2 quantities with the logarithmic de- 
pendence on Q 2 for fixed xp and (3. They are related 
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FIG. 7: A comparison of a® from the two considered dipole models with the newest ZEUS Collaboration data 20]. The 
solid lines correspond to the GBW parameterisation of the dipole cross section with the color factor modifications and 
(|12[) . while the dotted lines correspond to the CGC parameterisation. The dashed lines show the results without the charm 
contribution. 



to the diffractive parton distributions by the standard 
collincar factorisation formulae, e.g. in the leading log Q 2 
approximation we have: 



F 2 D ( X]P ,p,Q 2 ) 



E4 



(30) 



where <?j? and arc diffractive quark/antiquark distri- 
butions. We additionally assume flavour democracy for 
these distributions to account for vacuum quantum num- 



ber exchange responsible for diffraction, 



1 



2N 



f 



(31) 



where Y, D is a diffractive singlet quark distribution. This 
distribution is evolved in Q 2 by the DGLAP equations 
together with the gluon distribution g D . In contrast to 
the dipole model case, the xp dependence of the diffrac- 
tive parton distributions is fitted to data, as well as their 
form in (3 at some initial scale Qq. In Fig. [TO] we show 
the result of such an analysis (dashed lines) applied to 
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FIG. 8: The same as in Fig. [7] but for higher values of Q 2 . The dashed lines show the contribution without charm. 



the ZEUS data 0. In the small-/? region, the DGLAP 
fit curves are going through the experimental points with 
larger logarithmic slope, dF® /dlnQ 2 , than in the dipole 
approach. This illustrates the importance of more com- 
plicated diffractive states than the qqg state. 



VI. CONCLUSIONS 

We presented a comparison of the dipole model results 
on the diffractive structure functions with the HERA 
data. We considered two most popular parameterisations 
of the interaction between the diffractive system and the 
proton (the GBW and CGC parameterisations) which are 



based on the idea of parton saturation. The three compo- 
nent model with the qq and qqg diffractive states describe 
reasonable well the recent data. However, the region of 
small values of (3 needs some refinement by considering 
components with more gluons and qq pairs in the diffrac- 
tive state. This can be achieved in the DGLAP based ap- 
proach which sums partonic emissions in the diffractive 
state in the transverse momentum ordering approxima- 
tion. We extracted the diffractive gluon distribution from 
the dipole model formulae to use it for the computation 
of the charm contribution to . We found good agree- 
ment with the HERA data on the diffractive open charm 
production both for the the gluon distributions from the 
considered dipole models and the DGLAP fits from [28[. 
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FIG. 9: The same as in Fig. [7] but for the HI Collaboration data [l| . The dashed lines show the contribution without charm. 



The latter statement, however, is possible to make only 
due to present accuracy of the charm data. The results 
presented in this work might be a starting point for the 
future electron-proton collider LHeC at CERN. 



APPENDIX: DIFFRACTIVE GLUON 
DISTRIBUTION FROM DIPOLE MODELS 
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A comprehensive discussion of the derivation of the 
diffractive partem distributions in dipole models can be 
found in 31(. Here we only recall the derivation of the 
diffractive gluon distribution g D (xjp,z,Q 2 ) which sup- 
plements that in [3l|. We start from Eq. (jSJ) which 
we reduce to the collinear factorisation form. Let us 
substitute (1 — z)Q 2 — > Q 2 in there. We numerically 
checked that such a substitution practically does not 
change the diffractive structure function. Then the loga- 
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FIG. 10: A comparison of the GBW dipole model predictions (solid lines) and the results from the DGLAP fit (dashed lines) 
[H to the ZEUS data [|. 



rithmic derivative of F, 



reads 



On the other hand, from the DGLAP evolution equation 
we have for the diffractive singlet quark distribution (|3l"j) 



dR 



(qqg) 



81/3a s 



dlnQ 2 5127r 5 xjpB d 

z 



E4 f- 



{l-zf 



dk 2 (j) 2 2 



dY, 1 



a s {Q 2 ) f 1 dz 



d\nQ 2 



2tt 



x g D {x P ,z,Q 2 



N f 



I) 



(A.2) 



where we neglected on the right hand side a contribution 
(A 1) with the singlet quark distribution which is much smaller 
than the gluonic one when (5 <C 1. Thus from Eq. (150)) 
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FIG. 11: The diffractive gluon distribution xgo{xip,z,Q 2 ) 
for Q 2 = Am 2 c , 7.84, 10, 100 GeV 2 (from bottom to top) and 
xp — 0.0042. The solid lines correspond to the GBW param- 
eterisation while the dashed lines to the CGC parameterisa- 
tion. The dotted lines show the gluon distributions from the 
DGLAP fit [11] to the HI Collaboration data. 



wc find for the diffractive structure function 



dF 



D 




<91nQ 2 2tt - 

x g u (x P ,z,Q z ) 



(A.3) 



For small j3 we have: F^? q9 ^ « F®, thus by the com- 
parison with Eq. (|A.1[) . we find the following diffractive 
gluon distribution 



g D (x F ,z,Q 2 ) 



81 



256T: 4 x P B d (1 - zf Jo 



Q 2 



dk 2 <j> 2 2 



where 



(A.4) 



k 2 



dr r K2 (^\f zj (1 — z)krj Jz{kr) a(xp,r) . 



(A.5) 

For the GBW parameterisation of the dipolc cross sec- 
tion, wc additionally rescale the gluon distribution, 



1 



{Ca/CfY 



(A.6) 



and use formula (jlip for the dipole cross section. For the 
CGC parameterisation this rescaling has been already 
taken into account. In Fig.[Tl]wc show the gluon distribu- 
tions computed for the GBW parameterisation with the 
color factor modification (solid lines) and for the CGC 
parameterisation (dashed lines). There is practically no 
difference between them for the indicated scales. For 
Q 2 > Am 2 , the Q 2 dependence of the gluon distribution 
(|A.4|) is already very weak and close to the asymptotic 
limit obtained for Q 2 — > 00. We also show in this fig- 
ure the gluon distributions found in a DGLAP fit with 
higher twist to the recent HI data [28| (dotted lines) with 
a strong dependence on Q 2 due to the DGLAP evolution. 
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